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ABSTRACT 


The  integrals  which  determine  the  physical  parameters 
of  chief  interest  in  the  cavity  theory  of  Plesset  and  Shaffer 
are  systematically  considered.  Their  numerical  integrations 
are  shown  to  be  very  accurate  with  the  exception  of  the  values 
for  cavity  sine  at  low  cavitation  numbers.  Where  necessary, 
these  values  have  been  replaced  with  accurate  tabulations 
computed  by  use  of  suitable  series  expansions;  also,  asymp- 
totic formulas  valid  for  low  cavitation  numbers  are  presented. 


SYMBOLS 


a. 

b. 


C 


D' 


^3D* 


C 


P' 


D. 

d, 


d , 

m 


°n' 


'n* 


Al. 


J 

n 


I 


p. 

p. 

Pk' 


Q, 

R. 
s, 


length  of  wedge  face  (Fig.  1). 
parameter  fixing  scale  of  problem, 
two-dimensional  drag  coefficient: 

three-dimensional  drag  coefficient: 

Sd  = i P ''i) 

pressure  coefficient: 

^p=  <P‘Poo)/iP'^m  • 


drag  force. 

width  of  base  of  wedge  or  cone, 
maximum  width  of  cavity. 

Fourier  coefficient.  Sec  Eq.  (43). 
Fourier  coefficient.  See  Eq.  (44). 
Fourier  coefficient.  See  Eq.  (44). 
integral  defined  by  Eq.  (45). 
estimate  of  (I  , , - I ) when  n-*6. 
integral  defined  by  Eq.  (46). 
length  of  cavity  (Fig.  6). 


pressure  at  any  point, 
pressure  in  the  cavity, 
pressure  at  infinity. 


Q(^.  p)  = cos^  sin^  ^ 

integral  defined  by  Eq.  (49). 
integral  defined  by  Eq.  (50). 


t,  complex  variable.  Sec  Eq.  (2). 

v^,  velocity  at  infinity. 

v^,  ratio  of  velocity  at  infinity  to  that 
on  cavity  wall. 

W,  complex  potential;  W = U + iV. 

X,  abscissa  of  maximum  cavity  width 

(Fig.  1). 

y,  vertical  distance  from  A to  I 

(Fig.  1). 

dY.  dimensionless  coordinate  defined  by 

Eq.  (9). 

Y,  integral  defined  by  Eq.  (13). 

z,  complex  physical  coordinate: 

z = X + i y. 

a,  parameter  related  to  cavitation 

number  by  Eq.  (4). 

p,  wedge  or  cone  half-angle. 

A,  = a(;^,  q)  = (1  - sin^a  . 

A',  = A'(j(,  a)  ■ ( 1 - cos^a  cos*f()  ^ 

^ , dimensionless  complex  velocity, 
normalized  to  unity  on  AIA'. 

T] , complex  variable. 

8,  arg  (;  . 

X = 2 p/ir  . 

p,  fluid  density. 

<r,  cavitation  number. 

parameter  (real). 

argTj. 


Note:  The  notation  used  for  elliptic  integrals  E,  K,  etc.,  is  that  of  Ref.  8. 


IV 


INTRODUCTION 


The  usual  criticism  of  the  classical  wake  theory^  is  that  it  fails  to  account 
for  the  large  suction  pressure  on  the  downstream  side  of  the  body.  This  objec- 
tion is  in  part  overcome  by  the  theory  of  Riabouchinsky^,  who  constructed  a 
mathematical  model  of  the  wake  flow  around  a flat  plate  by  placing  an  image  plate 
downstream  with  a wake  region  at  constant  underpressure  in  between.  Whatever 
the  eventual  usefulness  of  his  theory  for  wake  flows  may  be,  it  is  now  certain 
that  such  a model  is  very  helpful  in  cavity  flow  analysis.  Here  the  underpressure 
in  the  cavity  is  specified  dimensionlessly  by  the  cavitation  number  c,  defined  by 


^oo  ^k 


(1) 


where  p and  v 

*^oo  GO 


are  the  pressure  and  velocity  at  infinity,  p^^  is  the  pressure 
in  the  cavity,  and  p is  the  liquid  density.  Riabouchinsky* s model  may  also  be 


Fig.  1.  The  cavity  flow  about  a symmetrical  wedge 
according  to  Plesset  and  Shaffer. 

useful  in  finding  two-dimensional  strut  shapes  which  have  small  tendency  to  cavi- 
tate, since  the  cavity  or  wake  may  equally  well  be  supposed  to  represent  a body 

with  constant  pressure  over  most  of  its  surface.  In  any  event,  Plesset  and 
3 4 5 

Shaffer  ’ ’ extended  Riabouchinsky' s theory  to  include  symmetrical  wedges 
(Fig.  1)  with  a view  to  applying  the  results  to  cavity  flow  problems.*  Their  work 


It  should  be  remarked  that  many  of  Plesset  and  Shaffer's  results  have  been 


given  in  an  interesting  report  by  J. 
declassified. 


W.  Fisher^,  which  has  only  recently  been 


I 


may  also  be  regarded  as  an  extension  of  the  results  of  Bobyleff  , who  considered 
the  symmetrical  wedge  with  an  infinite  wake,  which  corresponds  to  the  cavity 
flow  with  zero  cavitation  number. 

The  present  report  is  concerned  for  tne  most  part  with  the  evaluation  of 
certain  integrals  in  terms  of  which  the  cavity  parameters  of  interest  in  Plesset 
and  Shaffer's  theory  can  be  represented.  They  computed  these  integrals  numeri* 
cally,  in  most  cases  with  high  accuracy,  thereby  obtaining  reliable  values  for  the 
drag,  which  was  the  chief  concern  in  their  investigations.  The  cavity  dimensions 
at  high  cavitation  numbers  were  also  accurately  tabulated;  the  values  for  low 
cavitation  numbers,  however,  are  somewhat  in  error. 

Since  there  is  much  emphasis  at  present  on  applications  where  the  cavita- 
tion number  is  small.  Professor  Plesset  suggested  that  the  cavity  dimensions 
might  be  computed  analytically  in  this  special  case.  Such  calculations  have  in- 
deed proved  feasible  so  that  it  is  now  possible  to  evaluate  the  drag  and  cavity 
dimensions  with  reasonable  accuracy  for  any  case  likely  to  occur  in  practice. 

The  following  pages  give  a summary  of  these  calculations. 


THEORY  OF  PLESSET  AND  SHAFFER 


The  theory  of  Plesset  and  Shaffer,  which  uses  well-known  methods  of  con- 
formal mapping  to  solve  the  two-dimensional  potential  flow  problem,  will  be  out- 
lined here  for  reference  only.  For  a more  complete  exposition,  as  well  as  for 

remarks  on  the  comparison  with  experiment,  the  reader  should  consult  the  original 

,3,4,  5 
report  ’ ’ . 


In  Fig.  2 is  shown  the  hodograph-  or  ^ -plane*  which  can  be  constructed  by 
inspection  of  the  velocity  along  the  streamline  KCAIA'C'K'  in  the  physical-  or 
z-plane.  Also  indicated  is  the  plane  of  complex  potential  W = U f iV,  which  can 
be  drawn  immediately.  The  mathematical  problem  can  now  be  solved  by  con- 
formally mapping  the  ^ -plane  onto  the  W-plane.  The  mappings  which  effect 
this  are 


n 


t 


(2a) 

(2b) 


W 


b t 

(tan^o  + t*)^ 


(2c) 


* Plesset  and  t>.affer  normalize  the  velocity  to  unity  on  the  free  streamline 
AIA',  so  that  ^ is  ihr  ratfe  of  the  reflected  velocity  to  the  velocity  on  AlA*. 
Similarly,  the  potentia'  detined  by  Plesset  and  Shaffer  has  the  dimensions  of 
length,  rather  than  le^'g^ :..»)i^s  velocity. 
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where  t)  and  t are  auxiliary  variables  as  shown  in  Fig.  2.  These  mappings, 
along  with  the  fundamental  property  of  the  complex  potential. 


dW 

■ar 


t , 


enable  one  in  principle  at  least  to  calculate  the  potential  at  any  point  z. 


(3) 


t • » ♦ ly 

W U 

JL 

A ~~ 

c 

c 

2 - plane 


W « U+  iV 

1 

1 

1 

1 

1 

1 

K' 

C'  A' 

1 

• I 

A 

C 

K 

1 

1 

W- plane 


Fig.  2.  The  mappings  corresponding  to  Eq.  (2  a be). 


The  parameter  b fixes  the  scale  of  the  problem;  since  all  results  of  the 
theory  are  presented  in  dimensionless  form,  it  disappears  in  the  final  equations. 
The  parameter  o maybe  related  to  by  calculating  |^|  at  K in  terms  of 
W.  From  Eq.  (2a be)  one  finds  that 

V ^ 1 - sing 

o ” cos  a * 


But  application  of  Bernoulli's  equation  along  the  streamline  KCAIA*  in  con- 
junction with  Eq.  (1)  leads  lo 


and  hence 


1 + JT 


1 -f  sino 


rr 


1 + S’ 


cos  a 


(4) 
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Since  it  is  convenient  to  calculate  the  other  parameters  in  terms  of  a.  Table  1 
has  been  constructed  for  corresponding  values  of  3-  from  Eq.  (4). 

The  set  of  equations  (2abc)  and  (3)  enable  one  to  calculate  the  value  of  the 
velocity  at  any  point  in  the  s-plane,  and  the  shape  of  any  streamline,  including 
the  cavity  wall  AIA*.  Plesset  and  Shaffer  calculated  only  the  quantities  of 
chief  physical  interest,  namely,  the  pressure  along  CA,  from  which  the  drag 
may  be  obtained  by  integration,  and  also  the  coordinates  of  the  point  1,  defined 
by  X and  y*  (Fig.  1).  Each  of  these  calculations  will  now  be  considered  in 
turn. 


PRESSURE  DISTRIBUTION 

If  the  pressure  coefficient  is  defined,  as  is  customary,  by 


then  from  the  definitions  of  ^ and  v there  results 


(5) 


If  one  puts  t = sec  ^ along  the  wedge  face  AC,  the  pressure  coefficient  may  be 
expressed  as  a function  of  the  parameter  ^ by  use  of  Eq.  (2ab).  For  con- 
venience let  P)  = 1 - 1^1*  on  AC.  Then 


H4,  P)  = 1 


1 - ain^ 
cos  ^ 


so  that 

Cp  = (l+«r)  P(^p)  - (T  . 


(6) 


The  total  pressure  on  an  element  due  to  both  the  pressure  p on  the  front  and 
the  suction  pj^  on  the  back  is  given  by 

(7) 

In  order  to  find  the  point  to  which  the  pressure  p corresponds,  it  is 
necessary  to  integrate  dy  along  AC.  Now  from  the  relations  above. 


P-Pk 


= (l  + v)P(^p). 


d W .. 

^ ' TWr-TT 


$ 
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or  when  the  derivatives  are  evaluated,  one  has 


' "ICI  ti')^ 


dt  . 


(8) 


However,  along  AC  one  may  use  the  substitution  t = sec  ^ , Also,  in  what 
follows  it  is  convenient  to  consider  instead  of  d y the  dimensionless  variable 
dY  defined  by 


dY  = 

-(b  sinp  sin^Q  cos  a) 

dy  , 

(9) 

Then  in  terms 

of  ^ one  has 

dY  = 

P) 

a)  d^ 

(10) 

where 

P)  = 

cos^  sinf( 

1 -f  sin  ^ 
cos  ^ 

IF 

(11) 

and 

A(^.a)  = 

(1  - sin*a  sin*f()*^* 

• 

(12) 

Then,  if 

/a  — 

Y(^’.o.P)  = 

/ dY(^,o,p) 

(13) 

the  ordinate  of  the  point  in  question  is  located  by 


JUff 

sin  p 


(14) 


Therefore  by  picking  values  of  and  then  computing  and  y/a  sinp  from 
Eqs.  (6)  and  (14),  a plot  of  pressure  distribution  may  be  constructed.  An 
example  for  o’  = 0 is  shown  in  Fig.  3. 

While  Y(f^',  a,  p)  cannot  be  integrated  in  the  general  case,  several  special 
cases  can  be  evaluated.  For  p = ir/Z,  Pies  set  and  Shaffer  give 


Y(^,  a,  w/2) 


2 

sec  a 


cos  sinj^') 

(1  - sin^a  sin*^)'^* 


+ B(a)-  F(o,^)+ D(a,»(*) 


(15) 


where  B(a),  F(a,  ^),  and  D(o,  ^ are  elliptic  integrals  in  the  notation  of  Ref.  8. 
For  = 0,  this  reduces  to 


Y(0,o,w/2) 


1 

— \ 

cos  a am  a 


sin*a 


+ E - cos*Q  K 


(16) 
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Fig.  3.  The  pressure  distribution  on  a cavitating  wedge  (after 
Plesset  and  Shaffer).  The  plot  may  also  be  used  for  any  low  cavitation 
number  v by  considering  the  ordinate  to  be  a normalized  pressure 
coefficient  as  indicated. 


Also,  one  has 

Y(^',  0,  ir/2)  = cosf^'  + w/Z  - + sinj^'  cosf^'  . (17) 

!f  one  considers  the  case  a = 0,  p / v/2,  corresponding  to  Bobyleff's  problem, 
the  integrals  are  expressible  in  general  in  terms  of  the  logarithmic  derivative 

9 

of  the  gamma  fimction, 

Y(O.O.p)  = 4*1.  t 4! 

where 

^(*)  = log  T(z)  . 
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If  2P/t  is  a proper  fraction  2p/ir  * p/q»  say,  where  p and  q are  integers, 
the  integrations  may  be  performed  in  terms  of  elementary  functions  by  the  sub- 
stitution t = (1  • sin^)/(l  -f  sin^)  and  t = Then  in  terms  of  z,  the  integral 

is  a rational  function  which  may  be  integrated  by  standard  procedures. 

For  a>0,  the  binomial  expansion 

A = 1 + ”7  sin*a  sin^ff  + sin^a  sin^ ^ (19) 

may  be  employed  to  obtain 

Y(0,q,  P)  = Y(0,  0,  p) +-|8in*Q 

where  each  integral  in  the  series  is  now  integrable  by  the  same  procedure  as 
used  for  Y(0,  0,  p).  In  actual  practice,  it  appears  that  the  original  tabulations 
of  Plesset  and  Shaffer,  the  accuracy  of  which  can  readily  be  checked  by  compari- 
son with  these  special  cases,  will  prove  adequate  for  most  likely  applications. 
Their  results  are  shown  in  Table  2,  where  values  for  o = 1°  - 4°  have  been 
added  by  using  the  first  two  terms  in  Eq.  (20)  as  an  interpolation  formula. 


I 


Z 

sin*f(  0(4,  P)  df(  + 


(20) 


DRAG  IN  TWO  DIMENSIONS 

If  the  drag  coefficient  due  to  a drag  force  D is  defined  as 


(21) 


where  d is  the  width  of  the  wedge  base  (Fig.  1),  one  finds  by  integration  of 

the  total  pressure  on  each  element,  as  given  by  Eq.  (7),  that 

rr 


fjp(4.  P)  dY{4,  Q,  p) 
Cq  = (ff-+i)  V(6,q,P) 


(22) 


The  remarks  on  the  evaluation  of  Y(0,  a,  P)  above  also  apply  to  the  integral  in 
the  numerator.  In  general,  one  can  employ  the  expansion  Eq.  (19)  to  obt2un 
a series  of  known  integrals.  For  the  special  case  corresponding  to  Bobyleff's 
problem  the  integration  may  be  performed  exactly  to  obtain 

rf  2 

J^P(^.WdY{^.0.«  =4;^ 


(23) 


Apex  Half  Angl",  /9 


Fig.  4.  The  two-  and  three-dimensional  drag  coefficients 
at  e = 0 as  a function  of  p. 


and  for  Riabourhinsky's  case 


1 


T 

1 


P(^r/2)  dY{,^,o,ir/2)  = 


2 . 2 
cos  o sin  a 


E - cos'oK 


where  E and  K are  complete  elliptic  integrals.  Hence  for  these  cases  the 
drag  coefficients  are 


= 4 


A 


2 


9 


where  Y(0,  0,  p)  is  given  by  Eq.  (20)  and* 

Cjj(a,ir/2)  = 2(o-+l)  ^E(n)  - cos*a  K(n)j  . (26) 

As  before,  the  tabulations  of  Plesset  and  Shaffer  are  in  agreement  with  these 
formulas,  and  are  presented  here  for  convenience  in  Table  3 (see  also  Fig.  4). 


DRAG  IN  THREE  DIMENSIONS 


The  preceding  theory  is  not  applicable  to  the  three-dimensional  case  of  an 
axially  symmetrical  flow,  but  Plesset  and  Shaffer,  in  view  of  the  lack  of  any 
three-dimensional  results,  proposed  to  estimate  the  drag  for  cones  by  assuming 
the  distribution  of  pressure  from  center  to  edge  was  the  same  as  for  two  dimen- 
sions. It  is  now  certain  that  this  is  a reasonable  estimate,  useful  enough  for 
many  purposes,  although  a check  of  the  assumption  by  experiment  or  other  cal- 
culations would  be  of  interest. 

In  any  event,  when  the  pressure  is  integrated  on  the  cone,  one  obtains  for 
the  drag  coefficient  in  three  dimensions 


P(^p)  Y(^a,p)d  Y(^a,p) 

=o 

[Y(0,a,p)]^ 

The  integral  in  the  numerator  has  been  calculated  by  Plesset  and  Shaffer  for  the 
case  p = ir/2.  Using  Eq.  (15)  for  Y(j^,  a,  w/2),  they  obtain 


sec^o 


. 2 
sin  a 


■ t B'(a) 
2 sin  a tan‘a 


(28) 


In  the  event  a - a = 0,  the  integrations  may  be  performed  directly  or  the  limit 
of  the  above  expression  taken  to  find 


P(^,  p)  Y(^,  0.  n/2)  d Y{4,  0,  w/2) 


(29) 


so  the  zero -cavitation  number  drag  of  a disk  is 


* The  formula  for  given  by  Plesset  and  Shaffer,  Eq.  (25),  Ref.  4,  and 

Eq.  (25),  Ref.  5,  is  incorrect  by  a factor  of  2. 
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C3d(0.it/2)  - = 0.8053  (30) 

The  tabulations  of  Plesset  and  Shaffer  for  the  drag  coefficient  are  given  in  Table 
4 and  are  believed  accurate  enough  for  any  engineering  purpose  (see  also  Fig.  4). 


CAVITY  SHAPE  IN  TWO  DIMENSIONS 
The  location  of  the  maximum  width  of  the  cavity  (Fig.  1)  is  given  by 


where  0 = arg  ^ is  the  slope  of  the  cavity  wall  and  the  integration  is  to  be  per- 
formed along  the  free  streamline  from  A to  I.  Now  if  argt)  - ^ = vO/zp 
along  AIA*  one  has 


t ~ sin  ^ = sin(vO/Zp)  . 

Also,  on  this  contour,  since  | C | = If  then  |dz|  = |dw|  from  Eq.  (3)  and 
hence,  from  Eq.  (2  c)  and  the  above  value  for  t. 


where 


dz  = (cos a sin^a)bb*  cos^d^ 


A'  = (1  - cos*o  cos*  tjf  j 


2 j.  r*/2 


Then,  by  employing  Eq.  (9)  to  obtain  dimensionless  forms,  one  finds  that 


and 


X _ R(a,  P) 

a ■ ^(0,  o.P 


where  for  convenience  the  abbreviations 

w 

■7 


R(ttf  P)  = j'  cos  (>^  ^)  cos  A'  d«/^ 


(32) 


(33) 


(34) 
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and 


S{q,  p) 


z 

t 

sin 


tj/  ) sin  tf/  A'  dtf/ 


(35) 


have  been  used.  (If  one  wished  to  plot  the  shape  of  the  cavity  wall,  the  lower 
limits  would  be  changed  accordingly.  ) 

The  integrals  R(a,  p)  and  S(a,  P)  have  been  computed  by  expanding 
sin(Zpt)'/v)  and  cos(2p^/e)  in  Fourier  series,  and  then  integrating  term  by 
term.  Hence  the  construction  of  suitable  series  for  sin  and  cos  Xtf/  will 
now  be  considered,  where  X - 2p/tr. 


EXPANSION  IN  SERIES  OF  SIN  X^  AND  COS 

10 

If  X is  not  an  integer,  the  well-known  series 


VI  2X  . . 

cos  Xu/  = sin  X ir 

^ IT 


- + I (-1)"' 

2X  n>l 


cos 


is  convergent  and  may  be  integrated  term  by  term.  When  this  is  done  three 
times  (the  constants  of  integration  being  properly  evaluated),  the  following 
series  is  obtained: 

00 

sinX^/  = ^ X 

nmi  *' 

where  for  convenience  the  notation  has  been  used  that 


n«i  n (n  - X ) 


H.  = - — 
X IT 


2 V 

— sin  X IT  , 
w 


and 


D . = X*  (—  sinXx)  _LIL 

n»(n*-X*) 


(36a) 


(36b) 


(36c) 


It  is  now  necessary  to  find  trigonometric  series  for  ^ and  . These  can 
easily  be  constructed  from  the  series* 


12 


oo 


(2n+l)’ 


(37) 


n:| 


and 


= -12  f;  (-1)"^' 

n«l 


(38) 


since  inlegration  of  the  first  yields  an  expression  for  - 3irV^  » which  can 

then  be  solved  simultaneously  with  the  second  to  obtain  two  new  series,  one  for 
^ , and  one  for  . These  are 

oo 

(39) 

n»i 


and 


where 


00 


n<i 


B 


B 


^.11^  ± * -Sfl 


11,520  1 

rnmmmm 


7w^  n» 


$ n=  1 » 3|  5«  • 4 , 


H“2>  6*  • * 9 


(40) 


(41) 


and 


n 


n 


24x360  1 r,  , 2 . 37 X 

• l-(-l)  + “ -TTTr 

-2  9 ' ' nw  360 

7 X n L 


24x360 
7x*  * n“ 


, ■ n=l,  3,  5. . . , 


n— 2,  4,  6*  • • • 


(42) 


Then  the  series  for  sin  X can  now  be  written 


where 


oo 


sinXt//  = y g^sinn<^ 


ml 


(43) 


g =G.  B +H.  F + D.. 
®n  X n X n n X 


(43a) 


Also  the  above  series  for  sinX^  may  be  integrated  once  to  obtain 


13 


in  which 


and 


CO 


cos  X = h^  + X h^  cos  n 


n«i 

00 


8, 


h = 1 - X V _IL 
o ^_i  n 

n»i 


h = g Z*'  • 

n *n' 


The  values  for  and  h^  are  presented  in  Table  10. 


(44) 

(44a) 


(44b) 


EVALUATION  OF  INTEGRALS  R(q,  p)  AND  S(o,  p) 

The  problem  of  finding  the  integrals  R(a»  P)  and  S(a,  P)  can  now,  by  apply- 
ing Eqs.  (43)  and  (44),  be  reduced  to  summing  infinite  series  which  contain 
integrals  of  the  type 


and 


I 

n 


J 


n 


w 


T 

sinn 


f 

cos  n 


A'  ( ^ ) d 


A'(  V'  . 


(45) 


(46) 


By  writing  sinni^  and  cosn  tff  in  terms  of  powers  and  products  of  sin  tj/ 
and  cos  \j/  , the  integrals  I^  and  can  be  written  as  sums  of  integrals 

of  the  type 


/ 


sin**  ^ 


C0  8^  tff 


A'(  if/  ) d iff 


(47) 


Such  integrals  can  always  be  reduced  by  integration  by  parts  to  known  elementary 
functions  or  complete  elliptic  integrals.  For  the  present  work,  values  of  the 
integrals  I^  and  have  been  computed  for  a range  of  a,  and  the  results  are 
shown  in  Table  9.  Also  an  estimate  ot  Al  is  shown,  where 

.1  n>6  (48) 


since  this  quantity  is  used  in  a later  computation 
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With  the  values  of  A , b , I , and  J , it  is  easy  to  find  R and  S, 

n n n n 

since,  from  above,  there  results 


R(q,  p) 


p 

sin  a 


+ X 


09 

V h J 
n n 


(49) 


n»i 


and 


CO 


S(«.W  = Z 


(50) 


n«i 


These  series  were  computed  for  a sufficient  number  of  terms  to  assure 
reasonable  accuracy.  In  the  case  of  R(a,  p),  the  series  converges  very  rapidly 
so  that  only  five  terms  are  necessary  to  obtain  an  accuracy  of  ^ 0.  1%.  The 
series  for  S(a,  P),  however,  requires  at  least  seven  terms,  and  a further  correc- 
tion was  made  by  a procedure  as  follows. 

If  a is  small,  the  values  of  I„)  almost  equal  for  n >6  or  so, 

as  can  be  seen  from  Table  9.  Hence  for  an  improved  estimate  of  S we  write 

7 15  13 

S(o.p)  = 8n^n  + ^8  2 S (“-®)  «„  * 

nsl  n>«  nsS 

For  a ^ 5^,  the  resulting  values  are  accurate  to  within  4 0.  1%,  while  for 
a<  5°,  the  values  obtained  are  within  4 0.15%.  A summary  of  these  results 
is  given  in  Tables  5 and  6.  The  values  of  R(a,  P)  and  S(a,  p)  have  then  been 
used  with  those  for  Y(0,  a,  p)  from  Table  2,  to  find  x/a  and  y/a  ; the  results 
are  presented  in  Tables  7 and  8. 

For  the  case  of  p = w/2,  the  integral  R(a,  p)  can  be  expressed  in  terms 
of  elliptic  integrals  and  S{a,  p)  can  be  expressed  by  simple  functions,  as  was 
shown  by  Riabouchinsky.  There  results 


R(a.-;-)  = 


2 2 
sin  a cos  a 


E'(o)  - sin'o  K'(a) 


(52) 


and 


S(q,  = 


1 - sin  a 


cos  a sina 


(53) 


Then  for  Riabouchinsky' s case 


X _ E'  - sin^oK 

* sin^a  4 E - cos  aK 


(54) 
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and 


X 


sin  a 


1 - sing 

sin‘a  + E - coa®oK 


(55) 


The  values  from  these  equations  are  included  in  Tables  5,  6,  7 and  8. 

The  principal  cavity  dimensions  can  now  be  calculated  from  the  formulas 


and 


1 + c S c p X 
& 


""d" 


cot  P C S C P X 


(56) 


(57) 


These  parameters  are  plotted  in  Figs.  5 and  6.  In  Fig.  7,  the  ratio  g 
is  shown.  This  will  be  discussed  further  in  the  section  below  concerning 
asymptotic  values. 


APPROXIMATIONS  FOR  LOW  CAVITATION  NUMBERS 


If  the  cavitation  number  v is  small,  which  implies  that  a is-small, 
suitable  approximation  may  be  made  to  evaluate  the  several  integrals  discussed 
above.  For  example,  in  calculating  the  drag  coefficients,  one  may  use  the 
first  approximations 


and 


Y(0,o,p)  = Y(0,0,p)  + O(o^) 


Y(0,o,  p) 


E 

2 


PdY(0,  0,  p)  + 0(o“) 


(58) 

(59) 


so  that  the  drag  coefficient  becomes,  as  is  well  known, 

Cj3  = (l  + a)Cjj(0,p)  + 0(o*)  (60) 

The  same  approximation  holds  for  three  dimensions. 

For  computing  the  pressure  distribution  one  may  then  use  the  approxima- 
tion for  small  a 


(61) 


Dimensionless  Cavity  Length, 


.05  .10  .15  .20  .25 

P - Pk 

Cavitation  Number,  a = ^ / 

Fig.  7.  The  length>width  ratio  of  a two-dimensional  cavity. 

At  very  low  v,  the  ratio  is  independent  of  p,  being  asymptotic' to  Z/v. 
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If  instead  of  the  pressure  one  plots 


P-  Pl 


1+0-  i p 


= P(^.P) 


oo 


(62) 


a dimensionless  representation  is  obtained  which  is  exact  when  a ^ 0 and  is  a 
useable  approximation  for  any  small  cavitation  number  (Fig.  3).  The  accuracy 
of  the  approximation  may  be  noted  by  comparison  with  the  exact  values  in  Ref.  5. 
It  appears  that  the  error  will  be  less  than  a percent  or  so  if  a ^ 5°  (see  Table  1 
for  corresponding  values  of  <r  ). 


To  find  R(a,  p)  and  S(a.  P)  when  the  cavitation  number  is  small,  it  is 
only  necessary  to  obtain  approximations  for  sin  and  cos  Xtff  which  are  ac- 
curate near  the  origin.  For  this  purpose  one  may  write  sin  ^ X sin^f  . 
Then 


cos 


(63) 


so  that  for  a first  approximation 


sin‘^  ) cos  ^ ^'(  V'  • <^)  ^ V' 


= / cos  ^ h'(  ^ , a)  d - -J-  I sin*  ^ cos  }jf  A'(  ^ , o)  d ^ 

'^0  - 0 


1 


sin*  a 


1 + 


1 4 cos  a 


2 cos'^a  2 cos  a 1 - cos  a 


(64) 


For  o— *0,  we  have 


R(a.P)  = — T-  . 


sin*a 


(65) 


As  a-»0,  the  approximation  sinX^  = X sin  ^ may  also  be  used  to  evaluate 
S,  so 


S(< 


w 

7 


sin  iff  cos  if/  A' d ^ 


(66) 


I 


X 

~ sin  a 


i 
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Then  for  a— »0  the  cavity  length-width  ratio  becomes 


(67) 


But  for  o — • 0 we  have 


O'  * 2 X sin  a 


(68) 


so  that 


(69) 


Hence  at  low  cavitation  numbers  the  principal  cavity  dimensions  are  independent 
of  nose  shape,  being  functions  of  r only.  A similar  result  was  foind  by 
Reichardt^^  for  three-dimensional  cavities. 

It  will  be  found  that  the  approximations  given  above  will  enable  extension  of 
the  values  in  Tables  5,  6,  7,  and  8 to  lower  cavitation  numbers  with  reasonable 
accuracy. 
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TABLJC  1 - THT  CAVITATION  NUMBER.  *<a,  p) 


P«15° 

30° 

45° 

o 

o 

75° 

o 

o 

105° 

120° 

135° 

150° 

165° 

0» 

0.0000 

0.  0000 

0.0000 

0. 0000 

0.0000 

0.0000 

0.  0000 

0.  0000 

0.0000 

0.  0000 

0.0000 

t 

0.0058 

0.0II7 

0.0176 

0.0235 

0.0295 

0.  0355 

0.0416 

0.0476 

0.0538 

0.  0599 

0.0661 

2 

0.0117 

0.0235 

0.0355 

0.0476 

0.0598 

0.0722 

0.0848 

0.  0974 

0. 1103 

0. 1232 

0.  1364 

3 

0.0176 

0. 0356 

0.0538 

0.0724 

0.0913 

0.1105 

0.1301 

0. 1500 

0. 1702 

e.  1909 

0.  2118 

4 

0.0236 

0.«77 

0.0724 

0.0977 

0. 1235 

0.1500 

C.  1771 

0.2049 

0.  2333 

0.  2623 

0.2921 

5 

0.0296 

0.0600 

0.0913 

0. 1236 

0.  1568 

0.1910 

0.2262 

0. 2624 

0.2997 

0.  3381 

0.  3777 

10 

0.0602 

0. 1241 

0.  1918 

0.2635 

0.  3396 

0.4203 

0.5058 

0.  5965 

0.  69.'6 

0.  7945 

0.9026 

IS 

0.0923 

0.1931 

0.  3032 

0.4235 

C.5549 

0.6984 

0.8551 

1.0264 

1.2134 

1.4177 

1.6408 

20 

0.I26I 

0.2682 

0.4281 

0.6083 

0.8112 

1.0396 

1.2969 

1.5866 

i.9129 

2.2803 

2.6940 

2S 

0. 1622 

0.  3506 

0.5697 

0.8242 

1.  120i 

1.4639 

1.8635 

2.3279 

2.8676 

3.4948 

4.2237 

30 

0. 2009 

0.4423 

0.7321 

1.0801 

1.4981 

2.0000 

2.  6028 

3.  3267 

5. 1962 

3S 

0.2431 

0 . 54  5 3 

0.9210 

1.3880 

1.9685 

2.6902 

3.5873 

40 

0.  2896 

0.  6630 

1. 1445 

1.7655 

2.  5663 

3.  5989 

4S 

0.34IS 

0.  7996 

1.4142 

2.  2387 

3.3447 

tabu:  2 - THE  INTECRAl.  Y(0,a,  p) 


P 

a 

0-15° 

o 

o 

45° 

o 

0 

75° 

1 

o 

0 

105° 

120° 

135° 

150° 

165° 

0° 

0.  S9401 

0.71443 

0.87184 

1.0822V 

1.37195 

1.78543 

2.40471 

3.  39686 

5.15664 

8.87842 

20.55725 

1 

0.59415 

0,71461 

0.87208 

1.08260 

1.  37237 

1.78601 

2.40555 

3.  39809 

5.  15863 

3.88205 

20.56616 

2 

0.  59458 

0.71517 

0.87278 

1.08353 

1.  37362 

1.78774 

2.40801 

3.40180 

5.  16459 

8.89293 

20.59291 

3 

0.59531 

0.71609 

0.87396 

1.08509 

1.  37571 

1.79062 

2.41213 

3.40797 

5.  17452 

8.91108 

20.63746 

4 

0.59631 

0,  71738 

0.87562 

1.08727 

1.  37864 

1.  7946t 

2.4179) 

3.41660 

5.  18844 

8.93647 

20.69935 

5 

0.59760 

0.71904 

0.87774 

1.09007 

1.  38240 

1.79985 

2.42533 

3.42771 

5.  20632 

8.96913 

20.78006 

10 

0. 60858 

0.73300 

0.89582 

1.  11385 

1.41436 

1.84399 

2.48844 

3.52238 

5.  35898 

9.24832 

21.46678 

IS 

0. 62762 

0.75736 

0.92745 

1. 15575 

1.47115 

1.92315 

2.  60292 

3.69665 

5.64503 

9.78311 

22.81780 

20 

0.65529 

0.79276 

0.97331 

1.21622 

1. 55255 

2.03556 

2.76341 

3.93657 

6.  02969 

10.48055 

24.51212 

25 

0.69356 

0.84189 

1.03744 

1. 30136 

1.66807 

2. 19659 

2.99594 

4.28896 

6. 60392 

11.54  241 

27.15383 

30 

0.74451 

0.90768 

1.12366 

1.41646 

1.82512 

2.41681 

3.31584 

4. 77666 

7.40369 

13.03090 

30.88178 

35 

0.81149 

0.99468 

;.  23846 

1.57076 

2.03710 

2.71615 

3.75382 

5.44951 

8.51562 

15. 11691 

36.14913 

40 

0.89944 

1. 10981 

1.39153 

1.77819 

2.  324  53 

3. 12554 

4.  35817 

6.38651 

10.07890 

18.07856 

43.70293 

45 

i.OiSTi 

1.26347 

1.39791 

2.06056 

2.71977 

3.69456 

5.20743 

7.71805 

12.32599 

22.  38586 

54.82158 

tabu:  3 - DRAG  COETriCIENTS  FOR  TWO-DIMENSIONAL  WEDGES,  Cp(a.  P) 


« 

P«15° 

30° 

45° 

60° 

o 

1 

*40 

o 

0 

105° 

120° 

135° 

150° 

165° 

0° 

0.2838 

0.4885 

0.  6370 

0.7448 

0.8230 

0.8798 

0.9207 

0.9498 

0.9703 

0.9845 

0.9939 

5 

0.  2927 

0.5185 

0.6959 

0.8375 

0.9528 

1.0482 

1.  1293 

1.1993 

1.  2613 

1.3176 

1.3694 

10 

0.  3029 

0.5521 

0.  7624 

0.9443 

1. 1057 

1.  2522 

1.3885 

1.5180 

1.64  36 

1.  7674 

1.8914 

15 

0.3146 

0.5901 

0.8385 

1.0688 

1.2880 

1.5016 

1.  7142 

1.9297 

2. 1514 

2.  3825 

2.6259 

20 

0.3281 

0.6330 

0.9258 

1.2148 

1.5073 

1.8097 

2.  1279 

2.4676 

2.8346 

3.2346 

3.6741 

25 

0.3436 

0.  6826 

1.0278 

1.3889 

1.7754 

2.1966 

2.6621 

3.1825 

3.  7695 

4.4364 

5.  1975 

30 

0. 3618 

0.7401 

1.  1481 

1.5993 

2.  1080 

2.6898 

3.3634 

4.1497 

5.0738 

35 

0.  3829 

0.  8076 

1.  2922 

1.8574 

2.  5274 

3.3308 

4. 3031 

40 

0.4080 

0.8880 

1.4676 

2.  1803 

3.  0676 

4. 1830 

45 

0,437$ 

0,9$S3 

-- 

1.  6852 

2.5927 

3.  7806 

tabu:  4 - DRAG  COEFFICIENTS  FOR  CONES  OF  REVOLUTION,  Cq(«.P) 


• 

P«15° 

o 

o 

45° 

o 

o 

75'* 

90° 

105° 

120° 

135° 

150° 

165° 

0° 

0.2045 

0. 3758 

0.5181 

0.6350 

0.  7296 

0. 8053 

0.8646 

0.9101 

0.  9442 

0.9693 

0.9874 

5 

0.2108 

0. 3990 

0.  5664 

0.7143 

0.8450 

0.9599 

1.0609 

1.  1495 

1.2276 

1.2972 

1.  3604 

10 

0.21*4 

0.4255 

0.  6213 

0.8064 

0.9*17 

1.1477 

1.3054 

1.4559 

1,6004 

1.7407 

1.8792 

15 

0.  2275 

0.4557 

0.  6*49 

0.9147 

1.1459 

1.3787 

1.6140 

1.8527 

2.0966 

2.  3478 

2.  6097 

20 

0.2379 

0.4906 

0.  7586 

1.0429 

1.  3445 

1.6683 

2.0072 

2.  3725 

2.  7649 

3.1891 

3.6525 

25 

0.2501 

0.5311 

0.8456 

1.  1969 

1.5*92 

2.0274 

2.5171 

3.  0655 

3.6*17 

4.  3779 

5. 1686 

30 

0.2648 

0. 5788 

0.9494 

1.  3847 

1.8950 

2.4918 

3.  1898 

4.0061 

4.9635 

35 

0.2*18 

0.6355 

1.0749 

1.6173 

2.2*35 

3.09*8 

4.0950 

40 

0.  3022 

0.7039 

1.2296 

1.9110 

2.  7*82 

3.9109 

45 

0.  3270 

0.7175 

1.4236 

2.2*99 

3.4593 

TABLX  S - THE  INTEGRAi.  A(«.  Si 


P*15° 

30° 

GO 

OO 

3283.  1 

3283. 1 

821. 00 

820.89 

365.  06 

364.95 

205.48 

205. 37 

131.61 

131.52 

33. 139 

33.067 

14.909 

14.851 

a.  5323 

8.4838 

5.5853 

5.5447 

3.9877 

3.9512 

3282.9 

820.68 

364.76 

205.20 

131.  36 
32.948 
14. 795 

8.4038 


90° 

105° 

OO 

GO 

3281.2 

3280.  5 

819.41 

818.84 

363.68 

363.  17 

204.23 

203.79 

130.47 

130.07 

32.315 

32.021 

14.256 

14.027 

7. 9908 

7.8037 

5.1222 

4.9723 

3.5836 

).4460' 

1 

00 

3279.9 

818.22 

362.63 


135° 

150° 

165° 

GO 

GO 

OO 

3279.  5 

3279.  1 

3278.2 

817.63 

816.97 

816.14 

362.  05 

361.42 

360.  69 

202.  76 

202.  18 

201.52 

129.12 

128. 59 

128.00 

31.  337 

30.955 

30.551 

13.498 

13.  205 

12.898 

7.  3726 

7. 1352 

6.8870 

4.  6073 

4.4068 

4. 1988 

3.  1328 

2.9631 

2.  7885 

TABLX  6 • THE  INTECRAI.  S(a,  p) 


a 

P»15° 

o 

o 

45° 

0 

o 

75° 

0 

o 

105° 

120° 

135° 

isr° 

165° 

OO 

GO 

GO 

00 

OO 

00 

OO 

OO 

OO 

00 

GO 

1 

9.4323 

18.861 

28.  273 

37.  659 

47.01! 

56.  316 

65.  578 

74 . 790 

83.966 

93.  175 

102.  36 

2 

4.6529 

9. 3010 

13.936 

18.550 

23.  136 

27.  687 

3?.  204 

36. 678 

41.  116 

45.546 

49.939 

3 

3.0615 

6.  1179 

9.  1618 

12.187 

15.188 

18.  157 

21.  094 

23.994 

26.857 

29.  700 

32.  504 

4 

2.2668 

4.5235 

6.7784 

9.0108 

11.220 

13.401 

15.  551 

17.665 

19.744 

21.797 

23.811 

5 

1.7910 

3.5768 

5. 3514 

7.  1091 

8.8450 

10.  554 

12.  234 

13.879 

15.490 

17.073 

18.617 

1 

0 

0.84477 

1.6845 

2.5155 

3.3322 

4. 1306 

4.9067 

5.  6579 

6.3804 

7.0725 

7. 7340 

8.3608 

1 

5 

0.  53607 

1.0678 

1.5909 

2.  1013 

2.  5953 

3.  0693 

3.  5030 

3.9479 

4. 3480 

4.7200 

5. 0635 

2 

0 

0.  38504 

0. 76618 

1. 1396 

1.  5016 

1.8490 

2.1787 

2.4879 

2. 7749 

3.0384 

3.2772 

3.49»4 

2 

5 

0.29710 

0.59074 

0.87732 

1.  1535 

1.4164 

1.6633 

1.8918 

2.1002 

2.2874 

2.4524 

2.  5945 

3 

0 

0.  24052 

0.47777 

0.  70856 

0.92985 

1.  1389 

1.  3333 

1.5110 

1.6704 

1.8106 

1.9306 

2.0300 

TABLE  7 - THE  PARAMETER  s/«  WHICH  PUCES  THE  CAVITY  LENGTH  (Ft*.  1) 


220. 
54. 
23. 
20  I 1?. 
25  8. 


CO 

4594. 
1148. 
509.  6 
286.  3 

182.9 
45.  11 
.6 
. 7 
. 5 


GO 

3032. 
757.  1 


oo 

2392. 

597.0 

264.7 

148.4 


94.64 
23. 

9.1 
5. 

3. 159 


i 


90° 

0 

o 

GO 

00 

1837. 

1364. 

458.4 

340.0 

203. 1 

1 5S.  6 

113.8 

84.28 

72.49 

53.  63 

17.  52 

12.87 

7.  422 

5.  389 

3.926 

2.824 

2.  332 

1.660 

1.483 

1.039 

0.  983 

0.679 

0.668 

0.459 

965.2 
240.  5 
106.4 
59.50 


oo 

635.  7 
158.  3 
69.97 
39.08 


14.  34 
3.  347 
1.  350 
0.681 
0.  382 


TABLX  8 - THE  PARAMETER  y/a  WHICH  PUCES  THE  CAVITY  WIDTH  (Pig.  1) 


oo  { 00 

15.88  26.39 

7.826  13.01 

5.  143  8.544 

6.313 


00 

34.79 
17. 12 
11.23 

8 288 


1.715  I 
1.171  I 
0.846  I 


00 

27.26 
13.37 
8.  745 
6.431 


CO 

22.01 
10.78 
7.040 
5. 170 


CO 

16.28 
7.961 
5.  190 
3.805 


1. 

904 

0. 

836 

0. 

482 

0. 

313 

0. 

212 

0. 

M* 

TABLX  9 - THE  INTEGRA1.S  l^(a).  -Al(a).  AND  J^(a) 


• 

‘l 

*2 

‘3 

*4 

‘5 

*6 

‘7 

Al 

Jl 

^2 

^3 

56. 316 

111.  524 

165.084 

217.04 

267.  59 

316.80 

364.62 

46 

3281.  3 

3275.7 

3267.  5 

2 

27. 687 

54.  300 

79. 328 

102.84 

125.08 

146. 12 

165.95 

18 

819-42 

814.93 

808.42 

3 

18,  157 

35. 270 

50.859 

65.023 

78.020 

89.957 

100.73 

8 

363.  68 

359.  77 

354.27 

4 

13.401 

25.787 

36. 707 

46.283 

54. 796 

62.  367 

69.010 

5 

204.24 

200.  76 

195,  97 

5 

10.554 

20. 122 

28. 277 

35. 166 

4 1.084 

46.  169 

50.442 

2 

130.477 

127.  326 

123.024 

10 

4.9067 

8.9518 

11.  798 

13.726 

14.940 

15.  900 

16. 109 

-0.4 

32. 314 

30. 124 

27.  346 

IS 

3.0693 

5. 3772 

6.  6837 

7.2059 

7.4500 

7. 2652 

7. 1364 

-0.3 

14.256 

12.  580 

) 0.576 

20 

2. 1787 

3.6775 

4. 3150 

4. 3886 

4.2591 

4.0656 

3.7990 

-0.  3 

7.9910 

6.6307 

5. 1142 

25 

1.6633 

2.7136 

3.0134 

2.8717 

2.6356 

2.4177 

2.  1842 

-0.  2 

5.  1306 

3.9865 

2.8047 

30 

1.  3333 

2.  1088 

2.2222 

1.9862 

1.7284 

1.5364 

1.  3532 

-0.  2 

3.  5839 

2.6121 

1.6396 

35 

1. 1080 

1.7037 

1.7084 

1.4334 

1. 1868 

1.0326 

0.8971 

2.  6561 

1.8204 

0.8912 

40 

0.9470 

1.4196 

1.  3588 

1.0715 

0.8469 

0. 7289 

0.6308 

2.0859 

1.  3322 

0.6416 

45 

0.  8284 

1.2140 

1.  1127 

0.8260 

0. 6255 

0.5381 

0.4683 

1.6942 

1.0140 

0.41SU 

TABLX  10  - THE  COEFFICIENTS  g^(p)  AND  h^lp) 


n 

P-15® 

o 

0 

45® 

60® 

7S® 

105® 

120® 

135® 

150® 

165® 

1 

.282143 

.541363 

.756806' 

.911554 

.994212 

.931349 

.797843 

.615545 

.405599 

. 192216 

2 

-.084996 

-.  152410 

-.  186657 

-.175962 

-.113818 

. 158948 

. 34  98  1 6 

.553495 

.746517 

.903423 

3 

. 022288 

.039581 

.047670 

.043827 

.027361 

-.033885 

-,  067471 

-.091470 

-.095212 

-.067787 

4 

-.002656 

-.004763 

-.005831 

-.005488 

-.003532 

.004725 

.009783 

.013766 

. 014761 

.010691 

5 

-.000371 

-.000615 

-.000653 

-.000484 

-.000206 

-.000056 

-.000474 

-.001147 

-.001700 

-.001527 

6 

-.000350 

-.000627 

-.000763 

-.000723 

-.000465 

.000618 

.001272 

.001772 

. 001874 

.001332 

7 

.000471 

.000831 

.000991 

. 000900 

.000552 

-. 000651 

-.001249 

-.001621 

-.001596 

-.001059 

8 

-.000083 

-.000149 

-.000181 

-.000168 

-.000111 

.000146 

.000300 

. 000369 

.000438 

.000311 

9 

-.000065 

-.000112 

-.000127 

-.000110 

-.000059 

. 000049 

.000070 

.000055 

.000016 

-.000013 

10 

-.000027 

-.000049 

-.000060 

-. 000059 

-.000037 

. 000048 

.000099 

.000137 

.000147 

.000100 

11 

. 000097 

.000179 

.000225 

.000215 

.000141 

-.000185 

-.000373 

-.000497 

-.000498 

-.000331 

12 

-.000011 

-.000019 

-.000024 

-.000022 

-.000015 

. 000019 

.000038 

.000053 

.000057 

.000041 

13 

-.000018 

-.000030 

-.000038 

-.000033 

-.000018 

.000017 

.000027 

.000028 

. 000019 

.000008 

0 

.958946 

.840992 

.661122 

.442207 

.212099 

-. 167520 

-.270250 

-.297693 

-.251004 

-.  143633 

1 

.282143 

.541363 

. 7S68C6 

.911554 

.994212 

.931  349 

.797843 

.615545 

.405599 

. 192216 

2 

-.042498 

-.076205 

-.093328 

-.087981 

-.056909 

.079474 

. 174908 

. 276747 

. 373258 

.451711 

3 

. 007429 

.013194 

.015890 

.014609 

.009120 

-.011295 

-.022490 

-.030490 

-.031737 

-.022596 

4 

-.000664 

-.001191 

-.001458 

-.001372 

-. 000883 

.001181 

.002446 

.003441 

.003690 

.002673 

5 

-.000074 

-.000123 

-.000131 

-.000097 

-.000041 

-.000011 

-.000095 

-.000229 

. 000  340 

-."00305 

6 

-.000058 

-.000104 

-.000128 

-.000120 

-.000077 

.000103 

.000212 

.000295 

.000312 

.000222 

7 

.000067 

.000119 

.000142 

.000129 

. 000079 

-. 000093 

-.000178 

-.000231 

-.000228 

-.000151 

8 

••.000010 

-.000018 

-.000023 

-.000021 

-.000013 

.000016 

.000031 

.000039 

.000038 

. 000024 

9 

• , 000007 

-.000012 

-.000014 

-.000012 

-.000007 

. 000007 

.000011 

.000012 

. 000010 

.000006 

